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Abstract. During recent years some non-thermal radio emitting OB stars have been 
discovered to be binary, or multiple systems. The non-thermal emission is due to syn- 
chrotron radiation that is emitted by electrons accelerated up to high energies. The 
electron acceleration occurs at the strong shocks created by the collision of radiatively- 
driven winds. Here we summarize the available radio data and more recent observations 
for the binary Cyg OB2 No. 9. We also show a new emission model which is being 
developed to compare the theoretical total radio flux and the spectral index with the ob- 
served radio light curves. This comparison will be useful in order to solve fundamental 
questions, such as the determination of the stellar mass loss rates, which are perturbed 
by clumping. 



1. Introduction 

Several massive stars, of both Wolf- Rayet and OB-type, show evidence of non-thermal 
radio emission (lBiegingetalJll989h . This non-thermal radiation is produced by rel- 



ativistic electrons accelerated at collisionles s shock s through the Fermi mechanism. 
Following the theoretical model of Eichler & Usov ([1993^) the shocks are created by 



the strong radiatively-driven winds that collide in the binary or multiple syst ems. This 



model was confirmed by the observations of WR140 (iDougherty et al.l l2005). The ac- 
celerated electrons advect away cooling down due to adiabatic and inverse Compton 
losses. In the presence of the magnetic field, they emit synchrotron radiation in the 
radio band. The synchrotron radiation is partially absorbed by the free-free absorption 
mechanism, leading to flux variability locked to the orbital phase. This variability con- 
sists of two contributions: the intrinsic synchrotron variations due to the changes in 
collision energy linked to the orbital phase and the changing free-free absorption along 
the line of sight. 

Recently new observations for the 05 -i- 06-7 binary Cyg 0B2 No. 9 were ob- 
tained. Here these new radio data together with previously published ones are shown 
(Fig-[D- A new emission model, which is applicable both to Wolf-Rayet and OB stars, 
is presented. This model is being developed in order to compare the observations with 
the theoretical results. The comparison is essential to investigate important issues re- 
lated to massive stars astrophysics, such as the clumping effect on the mass loss rate 
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phase in 852.9-days cycle 



Figure 1. Cyg OB2 No. 9: 6 cm flux versus phase during the complete orbital 
period of 852 .9 days. The black points are the available VLA radio data from 
IVan Loo et al.l ( [2008 ), the red ones are our preliminary reductions of the recent ob- 
servations from EVLA. A sharp drop during periastron passage is clearly observable. 

determination of single stars, the particle acceleration mechanism at the shocks and the 
binary frequency of early-type stars. 



2. Observations: Cyg OB2 No. 9 



Van Loo et all ( 2008h discovered that the observed VLA radio fluxes of Cyg OB 2 No. 



9 at 3.6, 6, and 20 cm show recurrent variations with a period of w 2.355 years. They 
accordingly suggested that the star might be a binary system. Bi narity was establi shed 
soon afterwards by the discovery of line splitting in the o ptical (Naze et al. [2008h and 
a first orbital solution of Cyg OB 2 No. 9 was obtained by Naze et al. ' (•201C givmg a 
period of 852.9 days. The 6 cm radio fluxes from Van Loo et al. (2008) are shown in 
Fig. [T] as black points. The phase-locked flux variability, which is a fingerprint of non- 
thermal emission in a binary, is clearly observable. The measurements cover a period 
of about 20 years but with only a few data points around periastron passage. 

In the framework of a large multi-wavelength campaign initiated by one of us 
(YN) with the goal to monitor the 2011 periastron passage of Cyg 0B2 No. 9, the 
first observable since the discovery of binarity, we obtained new EVLA data. They are 
superimposed on the old ones as red points (the data reduction is only preliminary). 
The new data agree very well with the old ones and they show a very sharp drop just 
before periastron. 



3. Modelling Cyg OB2 No. 9 



A previous version of the emission code was developed and a pplied to the bina ry sys - 
tems Cyg 0B2 No. 8 A (iBlomme et allllOlOh and No. 9 (IVolpil(201 ll : [Volpi et alboTlh . 



Modelling of non-thermal radio emitters 
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The code presents several approximations: the input hydrodynamical variables were 
determined analytically, the orbital motion was not included, and the shocks were con- 
sidered coincident with the contact discontinuity. The strength of the shock was chosen 
to have a constant value. These approximations created a discrepancy between our 6 cm 
Ught curves, and the corresponding observations. In the case of Cyg 0B2 No. 9 our 
fluxes are too high and the maximum occurs too early. 

A new emission code is now being developed in order to solve the above discrep- 
ancy. The input hydrodynamical variables are obtained with the Athena code, evolving 
the gas hydrodynamics and including the orbital motion. 

The shock positions are found in the hydrodynamical results by looking for the 
density jump. The shock strengths are now determined point by point leading to a more 
consistent model for the particle acceleration. The electrons accelerated at the shocks 
are followed along the post-shock streamlines while they are moving away and cool- 
ing down due to the adiabatic and inverse Compton losses. The relativistic electrons 
in the presence of the local magnetic field emit radio synchrotron emission. The mag- 
netic field is locally estimated following Blomme et al. ( 2010.) . The total flux is then 
calculated including free-free absorption and emission. 

4. Future work 

The phase-locked variability in the total flux will be analyzed and compared with the old 
VLA and new EVLA data from Cyg 0B2 No. 9. The space of the model parameters 
will be better investigated in order to improve our physical understand of coUiding- 
wind binaries. The code will also be applied to other massive binaries and multiple 
systems. The complementarity between observations in different wavebands as optical 
and X-rays will provide relevant constraints on the modelling of Cyg 0B2 No. 9. This 
will help us to further understand the production of non-thermal radio emission in this 
system. 
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